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OFF-NUCLEAR AGN AS A SIGNATURE OF RECOILING MASSIVE BLACK HOLES 
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ABSTRACT 

During the final phases of inspiral, a massive black hole (MBH) binary experiences a recoil due 
to the asymmetric emission of gravitational waves. We use recent results from numerical relativity 
simulations together with models of the assembly and growth of MBHs in hierarchical cosmologies, to 
study the dynamics, statistics, and observability of recoling MBHs. We find that, at redshift z < 3, 
kicked non-rotating holes are typically found between 1 and 30 kpc from their galaxy centers, while 
rapidly rotating ones are typically between 10 and a few hundred kpc. A recoiling hole that carries 
an accretion disk may shine as an "off-nuclear AGN" while it moves away from the center of its host 
galaxy. We predict that, depending on the hole spin distribution and the duration of their active 
phase, a population of off-nuclear AGN may already be detectable at low and intermediate redshifts 
in present deep Hubble Space Telescope observations. The James Webb Space Telescope may discover 
tens of wandering AGN per square degree, most of them moving within their host halos on unbound 
trajectories. 

Subject headings: black hole physics - cosmology: theory - galaxies: nuclei - quasars: general 
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1. INTRODUCTION 

The massive black holes (MBHs) ob served today at 
the c enters of nearby galaxies (see, e.g.. iRichstone etaLl 
1998) are expected to have grown through a series of 
gas accretion episodes, when they shine as active galac- 
tic nuclei (AGN), and of binary coales cences following 
the merger of their host ga laxies (e.g., :Bcgelma n et al.l 
1980: Vol onteri et al.l l2003a) . Coalescing MBH pairs will 
give origin to the loudest gravitational wave events in 
the universe, and are one of the primary targets for the 
planned Laser Inte rferometer Space Antenna (LISA; e.g. 
ISesana et al.l [2004). In general, gravitational waves also 
remove net linear momentum from the binary and im- 
part a kick to the center of mass of the system. The 
outcome of this "gravitational rocket" has been the sub- 
ject of many recent numerical relativity studies. Non- 
spinning holes recoil with velocities below 200 kms -1 
that only depend on the binary mass ratio, wherear- 
eas much larger kicks are predicted for rapidly rotat- 
ing holes (ICampanelli et alj[2~007at iGonzalez et al.ll2007t 
iHerrmann et al.ll2007f r 

If it is not ejec t ed from the host altogether (e.g . 
Madau et all 12004 [Merritt et all [200l iHaimanl 12004 



Volonteri fc Reesl 120061; ISchnittman k. Buonannol l2007D . 



the recoiling MBH will travel some maximum distance 
and t hen return to the center subject to dynamical fric- 
tion (jMadau fc Quataert|[2004) . Historically, MBH ejec- 
tions h ave been proposed to explain the nature of ra- 
diojets (jSaslaw et al]ll974trValtonen fc Heinamakil2000l ) 
and t rails and fi l aments in spiral an d interacting galaxies 
(e.g. .lArd [19761 : lArp et al.ll200lD . lHaaue-Copilah et all 
(fi~99"7f T present a comprehensive compendium of early 
ideas on MBH ejections that have been recently revi- 
talized. 

Galaxy mergers are frequent at early times, so a sig- 
nificant number of coalescing MBH binaries is expected 
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to form then. Galaxy mergers are also a leading mecha- 
nism for supplying gas to their nuclear black holes, and 
a recoiling hole can retain the inner parts of its accretion 
disk, providing fuel for a continuing luminous phase along 
its trajectory. Two possible observational manifestations 
of gravitational-radiation ejection h ave then been sug- 
gested: (1) off- nu clear AGN activity (jMadau &: Quataertl 
l2(MlLo^l2lwl : and (2) broad emission lines that are 
substantially shifted in velocity relative t o the narrow- 
line gas left behind (|Bonning et al.ll2007t ). In this Let- 
ter we use models of the assembly and growth of MBHs 
in hierarchical cosmologies to study the statistics and 
dynamics of recoiling holes, and explore the conditions 
dictating their observability as off-nuclear AGN in deep 
optical and X-ray imaging studies. 

2. BASIC MODEL 

We follow the assembly history of dark matter ha- 
los and associated MBHs via cosmological Monte Carlo 
realizations of the merger hierarchy from early times 
to the present in a ACDM cosmology, using a model 
that captures many features of the MBH/ AGN popu- 
lation (e.g., luminosity function of quasars and its evo- 
lution, black hole mass density today, stellar core forma- 
tion due to MBH binary mergers). The main assump- 
tions of the models ha ve been discussed elsewhere (e.g. 
IVolonteri et"~aTI l2003allbl : IVolonteril 120071 and references 
therein). MBHs get incorporated through halo mergers 
into larger and larger structures, sink to the center ow- 
ing to dynamical friction against the dark matter back- 
ground, and form bound binaries. The MBHs in galax- 
ies undergoing a major merger (i.e. having a mass ratio 
> 1 : 10) have masses corresponding to the Mbh — re- 
lation of their hosts, i.e. at coalescence the MBH binary 
mass ratio scales with the central velocity dispersions of 
the progenitor halos. After the formation of the binary, 
the black hole pair shrinks and coalesces on a timescale 
that is shorter than the merger time of the progenitor 
halo s, as a result of stellar and gas-dynami cal processes 
(e.g. iBerczik et al.ll2006HMaver et al.ll2007l ). 

The recoil velocity Vjdck depends on the binary mass 
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ratio q = M2/M1, on the dimensionless spin vectors of 
the pair a\ and 0,2 (0 < a, < 1), and on the orbital 
paramete rs. All current num erical data on kicks can be 
fitted by (jBaker et al.ll2008f) 

Vkick = v m e x + v± (cos £ e x + sin £ e„) 4- U|| e z , (1) 

^^^^(l + B))), (2) 

U± =ffr / 2 (l + (Z )- 1 (4- (Za f), (3) 

V|| = Krf{l + g)- 1 cos(9 - 6 ) (at - qai) , (4) 

where 77 = q/(l + q) 2 is the symmetric mass ratio, the 
indices || and _L refer to projections parallel and per- 
pendicular to the orbital angular momentum, respec- 
tively, e x and e y are orthogonal unit vectors in the orbital 
plane, is the angle between (a 4 - — (Z^i ) and the sep- 
aration vector at coalescence, and Qq is some constant 
for a given mass ratio. Here, A = 1.35 x 10 4 kms _1 , 
B = -1.48, H = 7540 ± 160 k ms" 1 , £ = 215° ± 5°, 
and K = 2.4 ± 0.4 x lO^ms" 1 (jBaker et al.ll2008[ ). We 
will assume in the following that the orbital parameters 
and the spin orientations are isotropically distributed. 
The configuration producing the maximum recoil kick 
(equal-mass rapidly rotating holes with anti-aligned spins 
oriented parallel to the orbital plane, with the pair re- 
coiling along the z-axis ), Vkick = Krp = 3750 km s -1 , 
is then quite rare (c.f. iBogdanovic et afll2007f ). Little 
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Fig. 1. — Left panel: Kick velocity distribution for all z < 3 
MBH binaries with Mi + M2 > 10 J Mq in our cosmological real- 
izations. Solid curve: nonspinning holes. Dashed curve: rapidly 
spinning holes with a; = 0.9. The orbital parameters and spin 
vectors are assumed to be isotropically distributed. Right panel: 
distribution of binary mass ratios. Solid curve: nonspinning holes. 
Dashed curve: rapidly spinning holes. In the latter case more holes 
are ejected, and the distribution of binary mass ratios is slightly 
different. Dotted curve: Monte Carlo sampling the z = MBH 
mass function. 

is known about the masses of MBH binaries and their 
spins. The distribution of all binary mass ratios ex- 
pected in hierarchical models is found to be relatively 



flat (|Volonteri et al.ll2005D . In this work we focus on bina- 
ries that could give origin to detectable off-center AGN: 
we therefore select pairs coalescing at z < 3 with total 
masses Mi + M 2 > 10 5 M Q . The distribution of mass 
ratios for this subsample is shown in Figure [TJ For com- 
parison, we have also derived a distribution of mass ra- 
tios by Monte Carlo sampling th e z = MBH mass func - 
tion (adopting the Schecter fit in lMerloni fc H einz 2008), 
where we have assumed Mi,M 2 € [10 4 M Q , 10 10 M Q ], 
and we have imposed the condition Mi + M2 > 10 5 M . 
In the sampling of the z — mass function most binaries 
are composed of two rather small holes of similar mass, 
skewing the mass ratio towards unity ((q) = 0.17 ± 0.23, 
implying Vkick ~ 250 km s -1 for a.j = 0.9). The binary 
mass threshold we impose in the cosmological sampling 
also causes a loss of coalescing pairs with gwl, because 
of the limited major merger activity of the correspondin g 
host galaxies (of order of 0.1 Gvr" 1 . IWhite et al.|[2007t ). 

Large spins are a natur al consequence of prolonged disk 
accre tion episodes (e.g. iBardeenl Il970t iVolonteri et al.1 
2005), but can be avoided if AGN are fed by a se- 
ries of small-scale, ra ndomly oriented accretion events 
(|King fc P ringlc 2005|). To bracket the uncertainties, we 
consider here two extreme scenarios: prior to coalescence, 
both holes are assumed to be either non-rotating or are 
rapidly spinning with a = 0.9. The resulting distribu- 
tions of kick velocities for these two cases are plotted in 
Figure[TJ Close to 7% of all coalescing non-rotating holes 
are ejected from the nucleus with kick velocities between 
100 and 200 kms -1 ; in the spinning case, 12% get kicks 
above 100 kms" 1 and 2% above 500 kms" 1 . 

3. DYNAMICS OF RECOILING HOLES 

In our default model the host galaxy is described by 
two sp herical mass compone nts: a) the dark matter fol- 
lows a iNavarro et al.l (|1997h profile with median con- 
centration c med = 9(1 + z)- i (M halo /8 x 1O 12 /iM )-° 14 
and scatter around the median of Alogc = 0.14 
(Bull ock et al.l l2Q01h : 2) the central stellar profile is 
modeled instead as an isothermal sphere with one- 
dimensional velocity dispersion cr» and density 

2 

The velocity dispersion is related to the halo circular ve- 



locity V c at the virial radius by V c = \f 7 2,o* (cf. iFerrarese] 
2002). We truncate p* at the "bulge" radius Rb = 6kpc 
(cr*/200 kms -1 ) 2 : this ensures that the stellar bulge 
mass within Rb is equal to 1000 Mbh (jHaring fc Rixi 
l200l. using M RH = (10 8 M Q ) (<r*/200 kms" 1 ) 4 
(Trcmain e et al.l I2002T ) . The stellar profile is also trun- 
cated inside a core radius Rbh, as stars within the 
gravitational sphere of influence of the hole, Rbh = 
GMbh/c 2 , are bound to it and do not contribute to 
dynamical friction. The escape speed from such a 2- 
component potential can be approximated as w csc ~ 

2a[c 2 /(c)- 1 (l + c)- 1 + HR b /Rbh)] 1/2 , where /(c) = 
ln(l + c) — c/(l + c). Following a kick, the radial orbit of 
a MBH in a spherical potential is governed by: 



4ttG 2 In A p Mbh, 



-f(x)v (6) 



d 2 r _ GM (r) 

where f(x) = [erf(x) — (2x/^/Tr)e~ x ], x = v/y2a, and 
the velocity dispersion a is derived from the Jeans' equa- 
tion for the composite profile, assuming isotropy (e.g. 
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Fig. 2.— Left panels: Decay of a 10 6 M MBH in a spherical halo 
of mass Afhalo — 10"^ Mq and escape velocity v csc — 380 kms -1 , 
for VJtick = 0.8riesc (top) and = 0.6f e sc (bottom). Solid 

curves: orbits in a NFW+bulge potential. Dot-Dashed curves: or- 
bits in a NFW potential. Right panel: Distributions of MBH dis- 
placements for all z < 3 MBH binaries with M1 + M2 > 10 5 M Q in 
our cosmological realizations. Solid histogram: nonspinning holes 
in a NFW potential. Dotted histogram: nonspinning holes in a 
NFW+bulge potential. Dot-Dashed histogram: a = 0.9 holes in a 
NFW potential. Dashed histogram: a = 0.9 holes in a NFW+bulge 
potential. 



iBinney fc Trem ainc 1987|. Here M(r) describes the to- 
tal mass of the host galaxy within r, p(r) is the total 
density profile, and the second term represents dynamical 
fricti on against the stellar and dark matter background 
(e.g. iBinnev fe Tremaindll987l) . Stars and gas bound to 
the hole are displaced with it: here we assume, in first 
approximation (see below), that the total ejected mass 
is MsH.tot = 2Mbh- The Coulo mb logarithm, In A, in 
equat ion (|5J) is taken equal to 2.5 (|Gualandris fe Merritfj 
1*2007") . 

A MBH ejected with Vkick < v csc will oscillate about 
the nucleus losing orbital ene rgy via dynamical fric- 
tion (|Madau fc Quataertl 120041 ). Note, however, that 
the assumption that the central stellar density profile 
is unmodified by the heating effect of dynamical friction 
tends to artificially shorten the MBH decay timescale 
(|Bovlan-Kolchin et al.l20 04). Moreover, while in a spher- 
ical potential the hole will remain on a purely radial or- 
bit, in a realistic triaxial galaxy the hole acquires angular 
momentum and does not return through the dense cen- 
ter. This is known to reduce the effect of friction and to 
increase the de cay time by a factor of a few relative to the 
spherical case (jVicari et al J 1200*71 ) . To approximate the 
delaying effect of a triaxial geometry and of a decreasing 
core density, we have also run a case in which MBHs move 
in in a bulgeless pure NFW potential. The two left pan- 
els in Figure^ show two examples of orbital decays (for 
Vkick = 0.6,0.8w osc ) in a NFW+bulge ("short decay") 
and NFW only ("long decay") model. The right panel 
depicts the ensuing distributions of MBH displacements 
for all z < 3 kicked binaries in our cosmological realiza- 



tions, with and without spin. In all but the NFW+bulge 
a = case, a significant fraction of recoiling holes are 
found between f and 100 kpc from their galaxy centers 
(at redshift 2, a displacement of 10 kpc corresponds to 
an angular size of 1.2 arcseconds). Note that, in our re- 
alizations, only 0.1-0.3% of all MBHs at a given cosmic 
epoch are actually off-center. 

4. OFF-NUCLEAR AGN 

A recoiling hole that carries an accretion disk may 
be detected as it moves away from the center of its 
host. To assess the detectability of a subpopulation 
of off-nuclear AGN in deep imaging survey, we as- 
sume that all recoiling holes accrete at a fraction Je 
of the Eddington rate Me — ^GMsnm p /c(JT, L = 
cIeLe, with radiative efficiency e. The rest-frame op- 
tical part of their spectrum is modelled as a multicolor 
IShakura fc Sunvaevl (|1973l ) disk with maximum temper- 
ature fcT max ~ 1 keV (M BH / M©)" 1 / 4 , following a power- 
law with L v ~ i/ 1 / 3 at hv < fcTmax- The X-ray spectral 
energy distribution is described by a power-law with pho- 
ton index T = 1.9 and an exponential cutoff at 500 keV 
(Marconi et al. 2004). The duration of the luminous 
phase depends on the amount of disk material out to the 
radius R out ~ GMBn/V£ ick that is carri ed by the h ole. 
In the case of an a-disk, this is given by (lLoebll2007h 



M, 



disk 



(1.9 x 10 b M ) aZr(e^/f E )-^M^V 3 



(7) 



where e_i = e/0.1, M 7 = M B h/10 7 M , V 3 = 
Vkick/10 3 kms -1 , and a_i = a/0.1 is the viscosity pa- 
rameter. The condition < Mbh requires 



Viick > 550 kms" 1 aZt M {e-i/f E ) 



-0.21 M 0.43_ (8) 



When the above equation is not satisfied, we impose 
A-disk = A/bh corresponding to an AGN lifetime of 
«max = eca T /(4TtGm p f E ) « 4.5 x 10 7 yr (e_i// B ). 
A hole/disk system with (M7, a, e, /g, JVfdisk) — 
(1,0.1,0.1, 1, M BH ) recoiling with Vkick = 500 kms -1 , 
could then still be shining as far as 30 kpc away of the 
center of its host and appear as an off-center quasar. 

Figure [3] shows the expected number of detectable off- 
center AGN in deep and large surveys such as the Chan- 
dra Deep Field North and ffST-COSMOS, and the pre- 
dicted counts for future observations with the JWST. We 
set the minimum resolvable separation between a recoil- 
ing MBH and its galaxy centroid to twice the telescope 
resolution (i.e. 0.2 arcsec for ffSTand JWST, 1 arcsec for 
Chandra). We do not correct our counts for obscuration 
or Compton thick sources. In the HST-COSMOS fields 
(2 square degrees) we expect ~ 30 sources for the best 
case of large kicks (spinning holes), long decay timescales 
(no bulge), and long active phase (/b=Q.1, a=Q.01), but 
less than 1 in the unfavourable cases. In the small CDF- 
N field we expect at most 1 off-nuclear AGN in the most 
optimistic case, while in the 0.77deg 2 X-ray Chandra- 
COSMOS we can expect up to a few sources. Follow-up 
X-ray observations will be essential for the identification 
of the source as a genuine off-center AGN. The JWST 
yields the highest counts, about 40 sources per square 
degree in the best scenario. 

5. SUMMARY 
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Fig. 3. — Cumulative number counts of off-nuclear AGN de- 
tectable by the HST, Chandra, and the JWST. Solid curve: num- 
ber counts at the sensitivity limit of the HST-COSMOS survey. 
Dashed curve: same at the sensitivity of the CDF-N survey. Dot- 
ted curve: same at a it JWST sensitivity of 10 njy at 2fira. In 
each panel the upper (lower) set of curves shows the results for a 
bulgeless (NFW+bulge) potential. 

Motivated by recent numerical simulations of black 
hole binary coalescence and ensuing gravitational-wave 
kicks, and by the large rate of MBH binary formation ex- 
pected during the assembly of massive galaxies, we have 
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assessed the detectability of recoiling holes that retain 
the inner parts of their accretion disk along their trajec- 
tory and shine as off-nuclear AGN. While the search for 
kinematically offset QSOs in the Sloan Digital Sky Sur- 
vey (SDSS) indicates that large kicks rarely occur during 
a long-lasting active phase (with an upper limit on the 
incidence of wandering A GN with line-of-sight velocities 
above 800 kms" 1 of 0.2%. lBonning et al.ll2007t) . a strong 
candidate for a rapidly recoling MBH has been r ecentl y 
found in the SDSS database by iKomossa et all (|2008h . 
showing an emission line spectrum offset by 2650 kms -1 
(but see Bogdanovic et al. 2008, Dotti et al. 2008 for a 
different interpretation). Here, we have traced the for- 
mation history of MBH binaries in galaxies through cos- 
mological Monte Carlo realizations of the halo merger 
tree, assigned a kick velocity to every MBH binary coa- 
lescence expected in our simulation, followed the orbital 
evolution of the ejected hole in the potential of the host, 
and determined the duration of a off-nuclear AGN phase 
on scales that can be resolved by current and future op- 
tical, IR, and X-ray satellites. The basic conclusion of 
this Letter is that, depending on the hole spin distri- 
bution and the duration of their active phase, a popu- 
lation of off-nuclear AGN may already be detectable at 
low and intermediate redshifts in present deep HST ob- 
servations. The JWST should discover tens of wandering 
AGN per square degree, moving within their host halos 
on unbound trajectories. 
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